CMB Angular Power Spectrum from Correlated Primordial Fluctuation 
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The usual inflationary scenario predicts a Gaussian random primordial density fluctuation, differ- 
ent Fourier modes of which do not correlate with each other. In this paper we propose a correlation 
between these different modes. A simple case is that these different Fourier modes correlate with 
each other following a Gaussian function. For such a primordial density fluctuation we calculate the 
CMB angular power spectrum and flnd that its amplitude decreases but the decrease is different for 
different I. This feature can be used to constrain the the correlation strength from the real data. 

PACS numbers; 98.70.Vc, 98.80.-k 
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I. INTRODUCTION 

In the inflationary scenario, the primordial fluctuation 
comes from the quantum vacuum fluctuation of the infla- 
ton which drive the inflation 0] • Then the quan- 

tum fluctuation becomes classical when the correspond- 
ing scale goes out of the horizon. The primordial fluctua- 
tion produced in this method is regarded to be Gaussian 
random and scale invariant ^hich are confirmed by many 
astronomical observations [E B S [3 • 

This explanation of the origin of the primordial fluc- 
tuation seems reasonable and successful. However it is 
not a rigorous theory and needs more experimental tests. 
Besides the inflation, there exist other ideas to produce 
the primordial fluctuations, for example, the string gas 
model d, 113, [iH [3 ■ Even in the inflationary scenario, 
there are many kinds of modification to the primordial 
power spectrum, for example, the case that translational 
invariance is broken [l3| or rotational invariance is bro- 
ken [3, [H, m, [i3| • The quantum-classical transition in 
the inflationary scenario is also criticized by [l^. Con- 
sidering these uncertainties, our paper is motivated by 
checking one of the main features of the Gaussian random 
fluctuation from the experiments: the different modes of 
the primordial fiuctuation in the Fourier space are in- 
dependent from each other. As we know, the Cosmic 
Microwave Background (CMB) is the relics of early uni- 
verse, which can provide precise measurement on the 
primordial fiuctuation. So in this paper we introduce a 
small correlation between different modes of the primor- 
dial fiuctuation (in Fourier space) and then study how the 
correlation affects the CMB angular power spectrum. 

This paper is organized as follows. In section II we 
introduce the above mentioned correlations of the pri- 
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mordial fiuctuation in Fourier space. In section III we 
calculate the CMB angular power spectrum. The numer- 
ical results are presented in section IV and the discussions 
are given in section V. 



II. PRIMORDIAL FLUCTUATION 

The meaning of the Gaussian random primordial den- 
sity fluctuation is < >= P'{k)S^{k - k'), 
where ^(fc) is the Fourier mode of the primordial den- 
sity fluctuation £,{x). The delta function indicates that 
different mode does not correlate with each other. One 
can relax this property by allowing different modes to 
correlate. Generally speaking, one can rewrite the above 
equation as < ^{k)^*{k ) >= f{k,k ), where f{k,k ) is 
not required to be a delta function any more. How- 
ever f{k,k ) is not arbitrary, it must satisfy some con- 
straints since ^(x) is a real field. < £,{k)£^*{k ) >* = < 
£,{k')^*{k) > gives f*{k,k') = f{k',k), which means 
the real part of f{k,k ) is symmetric under the inter- 
change oi k,k , Ref{k, k ) = Ref{k , k), and the imag- 
inary part of f{k,k ) is antisymmetric under the inter- 
change of fc, fc , Imf{k, k ) = —Imf(k , k). Besides this, 
^*{k) = ^{~k), which indicates that < C{-k)(,{-k) 
f{-k',-k) = f{k,k') =< C(fc)r(^') >• Thus one 
obtains Ref{-k,-k') = Ref{-k',-k) = Ref(k,k') 
and Imf{—k,—k) = —Imf(—k,—k) = —Imf{k,k). 
P^{k)5^{k — k) is such a special function with a zero 
imaginary part. 

However, these constraints provide no positive infor- 
mation on the function f{k,k ). We must construct a 
concrete form of the above function. It is not strange 
to require that delta function is a good approxima- 
tion to this function, and noting that when cr — > 0, 

^ exp [— ^^2g"* ] ~* '^(2^~y)j 0116 can use a very simple 
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to replace the delta function. When ao — * 0, it ap- 
proaches 6^{k — k ). Here (Tq is a pure number and fco 
is the comoving wave number. This is a Gaussian form 
function, which is completely determined by the param- 
eter a = (TofcQ and its maximum is located at fc = fc . 
The FWHM is 2V21n2(T. Then in order to satisfy the 
above constraints, one can generalize P*(fc) to P^{V kk') 

or 

ps^k^y Their differ ence can be ignored in the weak 
correlation case (cr is small), but cannot be ignored in 
the strong correlation case (cr is great). However, it is 
easy to note that the strong correlation case is excluded 
by the CMB observation, so in this paper we only take 
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as an example to study the correlation effects on CMB 
angular power spectrum. In fact, the following discus- 
sions not only apply to primordial fluctuation with form 
(HI) but also to such a form /(fc, k ) = /(fc, k , \k — k |). It 
can be rewritten as 

f{k,k\\k-k\)^ I d'Kfik.k ,\K\)5^{k-k (3) 



III. CMB POWER SPECTRUM 

A. Two point correlator 

For scalar type of primordial density fluctuation, the 
CMB temperature anisotropy at direction n can be ex- 
pressed as 19, 20] 



T\n) 



(fkA^{riQ,k,k ■ h)S,{k) 



d^k I dr/5|,(77,fc)e-"^-"e(fc) (4) 





where S^{ri,k) is the transfer function, x = k{r]o — rf), 
•q is the conformal time and 770 is the conformal time of 
today. Following the standard methods, the temperature 
anisotropy can be expanded by the spherical harmonic 
coefficients 



dnYC^{n)T\n) 



(5) 



An important property of the modified primordial power 
spectrum < ^{k)^*{k') >— f{k,k',\k — k \) is that the 
CMB anisotropy produced by such a form of primor- 
dial density fluctuation is statistically rotational invari- 
ant which means that for any rotation R one always have 



(6) 



The definite proof of Eq.® is postponed to Appendix. 
With this property, one can obtain the following conclu- 
sion [U 

< a;m<V' >= ^ll'^mrn'Cl- (7) 

So the angular power spectrum C; is sufficient to reflect 
all the two point correlators of CMB anisotropy. 



B. Power spectrum 
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The CMB angular power spectrum can be calculated 



C, 
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( J dnYC^{n)T{n)) > (8) 

Substituting Eq.(I3|) into Eq.®, one obtains 
(4^)2 



Ci - 
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d^kd^kGi{k')Gi{k) 



d^AFiA), (9) 



has been used and F{A) is defined by 



where the formula e^*^'^ = ^T^J2i„i'i-''ji{kx)Yi*^{k)Yi„i{x) 



F(A) = I d^kd^kGiik')Gi{k)YUk') 

rn 

Y;^{k)f{k,k',\A\)S'{k-k - A), (10) 

where Gi{k) = dr]S^{ri,k)ji{x). In the Appendix we 
have proved 

F{RA)^F{A). (11) 
So, F{A) is only a function of A ee |A|. Thus 

Q =47r / A^dAF{A). (12) 
Jo 

Introducing new variables p — k + k ,p — k — k , and 
F(A) = 



one obtains 



8{2l + 



TtE I d^pd''pGi{k)Gi{k)Yi„,{k) 



YC„,{k)f{k,k,\A\)5'{f ~A) 
= Sfiri)^ / d'pGiik')Giik)YUk ) 

YC„,Ck)f{k,k,A) 
= \ f d^pGi{k')Gi{k)Pi{k' ■k)f{k,k',A),{n) 



3 



P + P 



P-P 



where k = — ^ — ,k = — ^ — . Since -F(A) is only 
a function of A, one can safely set A in the z di- 

^pI+pI + {pz + ^Y 



rection A 
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Noting that the integrand is only a function of + Pyj 
and using 
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one obtains 
Ci = 

Pi{k' ■k)f{k,k,A) 



dp, / pdpGi{k)Gi{k) 

-00 Jo 



(15) 



Finally noting that the integrand is an even function of 
Pz, one obtains a simple form, 
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and 



k' ■k = 
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IV. NUMERICAL RESULTS 

In this paper we use the Boltzmann code CMB- 
FAST ([l^l) to do numerical calculation and adopt a 
flat AGDM model with cosmological parameters: J7f, = 
0.045, Qc = 0.225, Qa = 0.73, Hq = 70, r = 0.08, = 1 
and A| = 1. Where t is the optical depth, n^. is spec- 
tral index and A^ is the primordial scalar perturbation 
amplitude. In CMB power spectrum calculation A^ is 
only a normalization factor, the value of which is deter- 
mined by observational data. Thus we can safely set it 
to be unit (usually the value of A^ is taken to be about 
2.35 X 10~^°), which will not affect the following discus- 
sion, since we mainly care about the relative change of 
the power spectrum. 

First we calculate {Gi{k)f = {J^^° drjS^ir], k)ji{x)f 
for I = 5,100,1000,2000 as a function of k. 
When there is no correlation between the different 
modes, the CMB angular power spectrum is Gi = 
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FIG. 1: Dependence of {Gi{k))^ on k, from left to right 
these curves correspond to (G5(fe))^34 * (Gioo(fc))2, 5500 * 
(Giooo(fc))2,37000* (G2ooo(fe))'' respectively. 
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FIG. 2: Dependence of G5(fc) on k. 



{Ai^f k'^dkP'{k){Gi{k)f, which reduces to Gi = 
(47r)^A| Jp°° dhik{Gi{k))'^ when the spectrum is scale in- 
variant. From Fig.([T|) it is easy to note that for different 
I the main contribution to Gi comes from different region 
of k. The larger is the main contributions to Gi comes 
from larger k. Then we plot G/(fc) ~ /J'" drjS^irj, k)ji{x) 
for I = 5, 1000, 2000 in their main region (see Fig.lgEED). 
One note that when I = 5 Gi{k) is almost positive except 
for small k where there exist a negative peak. One also 
note that for / = 2000 the average period is longer than 
that for I = 1000. 

The CMB angular power spectra are plotted in Fig. ([5]) 
and the relative angular power spectra changes and the 
relative cosmic variance are plotted in Fig.([n]). It is easy 
to note that the amplitude of angular power spectrum 
decreases when correlation is introduced. The greater 
is CT, the lower is the amplitude. From Fig.® it can 
be seen that the angular power spectrum decreases more 
significantly at the neighborhood of / = 1000 than at 
smaller or larger /. This is because of the phenomenon we 
mentioned in the last paragraph. Numerical calculation 
tells us that the contribution of F{A) to Gi (Eq.jn])) 
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FIG. 3: Dependence of Gioaa(fc) on k. 



0.00015 
0.0001 
5e-05 


-5ti-05 
-0.0001 
-0.00016 



I i I . 



0.14 0.145 0.15 0.155 0.16 0.165 0.17 0.176 



FIG. 4: Dependence of G2ooo(fc) on k. 



concentrate at the neighbourhood of A = v2(t and 
k' ■ k ~ 1 \ii this case when a is small (so Pi{k' • fc) w 1). 
Thus one can simply consider such a quantity Ci = 
(47r)2 dlnkGiik) * Gi{k + 5k) (refer to Eq.dHD) and 
compare it with Ci = {4:Tt)'^ d\nk{Gi{k))^ to under- 
stand the above mentioned phenomenon, where dk is a 
small fixed value. From Fig. (|2l3l4p one know that Gi{k) 
oscillate with fc, so one can further consider such a func- 
tion, 




X < -1 
-1 < X < 
< X < 1 
X > 1. 



The quantity Di = J f^{x)dx = |, while another quan- 
tity D2iS) EE J fix) fix + 6)dx = I - (5^(1 - I) (here 
5 <C 1). One can see that D2 < Di. As the same reason, 
the angular power spectrum decrease when the primor- 
dial density fluctuation is correlated in Fourier space. At 
the same time, when < Si < 62, D2iSi) > D2iS2)- Due 
to this, for a fix 6k, Ci decrease less in the case that Gi (fc) 
oscillates with a longer average period than the decrease 
in the case that G;(fc) oscillates with a shorter average 
period. Thus the angular power spectrum at large / de- 



FIG. 5: Unnormalized CMB angular power spectrum Ci as 
a function of I. The solid line corresponds to the Gaussian 
random case, while the long dashed line, the short dashed line 
and the dotted line corresponds to cr = 10"^", 5 * 10"^" and 
10~® respectively. 
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FIG. 6: Dependence of on I. 



ACi is defined as the CMB 



angular power spectrum with correlation minus the one with- 
out correlation. The solid line, long dashed line and short 
dashed line correspond to a = 10"^", 5 * 10"^" and 10~^ 
respectively. The dotted line is the c osmic varia nce of the 
angular power spectrum ACj/C; = i/2/(2? + 1) in the un- 
correlated case (here — ^2/(2i -I- 1) is plotted). 



crease less than that at about I ~ 1000, since Giooo(^) 
oscillates with a shorter average period compared with 
G20Qoik). At last one note that when Gi(fc) is always 
positive, the integrand Gi (fc) *Giik + 6k) is also positive, 
while when G;(fc) is positive and negative alternate, the 
integrand Gi(fc) * Giik + 6k) is not always positive. So 
when G;(fc) is always positive G; tends to decrease less 
than the case that G;(fc) takes positive value and nega- 
tive value alternately. This is the reason that the angu- 
lar power spectrum decrease less at small / than that at 
I = 1000. These features can be used to constrain the 
correlation strength parameter a from the real data. 
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V. DISCUSSION 

In this paper we only discussed a special type of the 
correlation of the primordial density fluctuation in the 
Fourier space. There exist other forms of correlation. 
From the point of experimental view, the correlation dis- 
cussed in our paper is typical because the correlation be- 
tween two modes is large only when they are almost the 
same. However, it is possible that the correlation in the 
Fourier space is scale dependent. Such possibilities can 
not be excluded. A better method is to divide these pos- 
sible forms into several types, and then consider their 
effects on CMB. Besides CMB temperature, the effects 
on CMB polarization and large scale structure also need 
a further study. 



Replace the dummy variables k, k' by i?fc, Rk' and note 
that S{Rk)^^ (fk, \Rk\ =}k\, d^{Rk') = d^k' , |i?fc'| = 
\k I and \R{k — k )\ = \k — k |, one obtains 



< T'{Rn )T'{Rn) > 

d^k d^k dT] S^{T],k)e" '^'^ 



dijS^iji, k)e 



-ixRk- Rfi 



f{k,k ,\k-k\).iAS) 



Since Rk ■ Rn — k ■ fi and Rk ■ Rfi = k ■ h, one finally 
obtain 
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< T'{Rn )T'{Rn) >=< T'in )T'{n) > . (A4) 



APPENDIX A: ROTATIONAL INVARIANCE 

In this Appendix we shall give a proof for the statisti- 
cally rotational invariance oiT^{n )T^{n) and rotational 
invariance of -F'(A). 

From Eq.Q, < T^{n )T^{n) > can be expressed as 



d^k J d'^k J dr]' Sl^{r],k )e 
d7jSl^iT],k)e''''''-^ <^{k)C{k ) > 



»)0 



\^ix k -h 







d^k / d^fc / dr]' S^{r],k)e 



'\^ix k -n 



no 



drjS^i?], fc)e"*=''=-"/(fc, k ,\k-k\). (Al) 



Then the two point ensemble statistical average after ro- 
tation R is 



The same reasoning can be applied to -F'(A). From 
Eg. (fTOj) . one obtains 



F{RA) = I d'kd'kGiik')Giik)YUk) 

m 

YC^{k)f{k, k,\A\)5^{k-k - RA). (A5) 



Replace the dummy variables k,k' by Rk, Rk' as 

before and noted that d^{Rk) — d^k, \Rk\ = 

|fc|, d^{Rk') ^ d^k', \Rk'\ = \k'\, \RA\ ^ 

|A|, Yir,.{Rk') = Ern'D'^^'imrnik'), Y*jRk) = 

J:m'D'*,{R)YC^{k) and S^[R{k - k' - A)] = 



\R\ 



6 {k ~ k — A), one gets 



< T'iRri )T''{RTi) > 



d^k / d^k / dr] S^{r],k)e'^ 



dr]S^{r],k)e''^'-''''f{k,k',\k-k'\). (A2) 



F{RA) = F{A). 



(A6) 
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